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ABSTRACT 


The  heat  transfer  associated  with  unsteady  film  con¬ 
densation  is  analysed  for  a  vertical  plate.  The  unsteady 
state  is  considered  to  be  created  by  the  time -dependent 
variation  of  either  the  uniform  wall  temperature  or  the 
g -force  field.  From  the  perturbation-type  theoretical 
treatment)  a  set  of  universal  variables  and  functionals 
derived  which  describes  the  unsteady  behavior  of  the 
film.  The  universal  functions  are  evaluated  and  tabu¬ 
lated  so  that  the  deviation  of  the  heat  transfer  from  the 
instantaneous  steady  state  value  can  be  computed  read¬ 
ily.  The  effect  of  varying  the  fluid  properties  on  the 
unsteady  heat  transfer  is  discussed. 
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SYMBOLS* 


g 

H 

h° 

k 

n 

Pr 

q 

T 

T 

V 

t 

u 

V 

X 


y 


/  lA 

(g/4|/  ) 


Specific  heat  at  constant  pressure 

Specific  heat  of  vapor  at  constant  pressure 

Steady-state  nt^ndimensional  stream  function 

Unsteady-state  nondimensionad  stream  function  defined  by 
Equation  (14) 

Perturbed  nondimensional  stream  functions 

Acceleration  body  force 

Steady-state  nondimensional  temperature 

Heat  of  condensation  of  vapor 

Thermal  conductivity 

Positive  integer 

Prandtl  number 

Local  heat  transfer  to  wall 

Absolute  temperature 

Absolute  temperature  of  vapor 

Time 

The  x-component  of  velocity 
The  y -component  of  velocity 

Direction  and  distance  along  plate  measured  from  leading 
edge 

Direction  and  distance  normal  to  plate  measured  from  wall 


Note: 


All  symbols  are  for  liquid  film  unless  specified  otherwise. 


iv 


r 


8 

e 


V 


Ifn) 


P 

0 

Superscript 

( 

Subscript 

s 

St 

w 


Steady-state  nondimensional  film  diickness 

Unsteady-state  nondimensional  film  thickness  defined  by 
Equation  (12) 

Perturbed  nondimensional  film  thickness 
Film  thickness 

Nondimensional  similarity  variable  defined  by  Equation  (10) 

Unsteady-state  nondimensional  temperature  defined  by 
Equation  (15) 

Perturbed  nondimensional  temperature 
Kinematic  viscosity 

Infinite  set  of  nondimensional  variables  defined  by 
Equations  (19)  and  (20) 

Density 

Stream  function  defined  by  Equations  (13) 

Total  differentiation  with  respect  to  the  variable  concerned 

Liquid-vapor  interface 
Instantaneous  steady  state  value 
Wall 


v 


INTRODUCTION 


Film  condensation  is  one  of  the  basic  heat  transfer  processes.  The  steady 
state  problem  has  been  widely  studied  since  the  time  of  Nusselt. 

Sparrow  and  Siegel^ analysed  the  relaxation  period  of  the  condensation 
process  following  a  sudden  drop  of  the  wall  temperature  below  the  condensation 
temperature.  The  analysis  was  based  on  an  approximate  method  which  neglects 
the  convection  terms  In  the  governing  eouatlons  and  assumes  a  linear  tempera¬ 
ture  profile  across  the  liquid  film. 

Oftent  In  engineering  applications!  the' condensation  takes  place  under 
conditions  which  are  continuously  unsteady.  The  unsteady  situation  is  usually 
caused  by  a  time -dependent  variation  of  the  wall  temperature.  With  the  current 
interest  in  space  technology!  a  factor  in  addition  to  unsteady  wall  temperature 
may  become  important  to  unsteady  condensation --that  is  the  unsteady  force 
field  which  prevails  in  a  space  vehicle  during  acceleration  and  deceleration! 
and  which  also  prevails  during  reduced-g  experiments. 

In  the  present  paper!  the  complete  boundary  layer  equations  will  be  used 
to  study  the  unsteady  film  condensation  process  and  accompanying  heat  transfer 
for  a  plate  located  parallel  to  the  acceleration  field.  The  analysis  applies  when 
either  the  uniform  wall  temperature  or  the  acceleration  field  is  arbitrarily  time 
dependent. 


Superscript  numbers  in  parentheses  refer  to  similarly  numbered 
references  in  bibliography  at  end  of  paper. 
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FORMULATION  OF  PROBLEM 


The  physical  model  studied  is  shown  on  Fig.  1.  The  plate  is  in  contact 
with  uniform  saturated  vapor*,  and  the  plate  temperature  is  considered  to  be 
always  below  the  saturation  temperature.  The  vapor*  therefore*  continuously 
condenses  at  the  liquid-vapor  interface.  The  heat  of  condensation  is  transferred 
to  the  plate  across  the  condensate  film  which  flows  along  the  plate  under  the 
influence  of  the  body  force.  The  process  is  unsteady  here  because  either  the  wall 
temperature  T^  or  the  body  force  g  is  arbitrarily  time  dependent.  In  the  present 
paper*  we  are  interested  in  studying  the  unsteady  behavior  of  the  liquid  film  so 
that  the  heat  transfer  to  the  plate  may  be  obtained. 

It  is  assumed*  in  the  usual  manner*  that  the  condensation  process  is  con¬ 
trolled  by  the  flow  within  the  liquid  layer  and  is  not  limited  by  the  supply  of 
vapor  at  the  interface.  That  is*  the  vapor  reacts  instantaneously  to  the  inflow 
requirements  defined  by  the  liquid  flow  conditions  at  the  liquid-vapor  interface. 

The  following  boundary-layer-type  conservation  equations  describe  the 
behavior  of  the  liquid  film.  The  density*  viscosity*  and  the  Prandtl  number  of 
the  liquid  are  assumed  to  be  constant  and  only  the  case  of  laminar  liquid  flow  is 
considered. 

Continuity  equation: 


du  8v  _ 
Sjc  9y 


0 


Momentum  equation: 


du 

Ft 


8u 


+ 


(1) 


(2) 


that  c  (T 

pv  V 


The  analysis  is  also  applicable  when  the  vapor  is  superheated  provided 

-  Tg)  «  h°. 
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Energy  equation: 


»T+„|T  . 

9 1  8x  8  y 

1/  8  ^T 

Pr  2 

d  y 

(3) 

The  boundary  conditions  are  as  follows: 

For  y  =  0 

u  af  0 

(4) 

V  >  0 

(5) 

T  «  T^(t) 

(6) 

for  y  =  5(xi  t) 

0 

ay 

(7) 

T  *  T^ 

(8) 

BY 

(9) 

Boundary  condition  (7)  states  that  the  shear  at  the  liquid-gas  interface  is 

negligible.  A  recent  study  of  steady-state  condensation  (3)  which  included  the 

Interface  shear  showed  that  the  boundary  condition  (7)  is  sufficiently  accurate  for 

most  cases:  the  only  exception  being  the  case  where  the  Prandtl  niunber  is 

c  (T  -  T  ) 

extremely  small  and|  at  the  same  time*  the  parameter  P  *  ^  large. 

Boundary  condition  (9)  is  derived  from  the  consideration  that  the  heat 
released  at  the  Interface  by  the  condensing  vapor  is  carried  away  into  tiie  film  by 
conduction.  The  momentum  equation  is  coupled  to  the  energy  equation  only 
through  the  boundary  condition  (9). 
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SOLUTION  OF  EQUATIONS 


In  the  present  sectiont  we  shall  solve  the  governing  equations  for  the 
velocity  and  temperature  profiles  through  the  film,  and  for  the  film  thickness  6  . 
The  results  will  then  be  used  to  obtain  the  heat  transfer. 

(21 

A  study  of  the  transient  analysis  given  by  Sparrow  and  Sieger  '  suggests 

that  the  response  of  the  film  to  the  time -dependent  variation  of  the  boundary 

condition  may  be  very  fast.  Therefore!  we  shall  obtain  an  unsteady  solution  In 

the  form  of  perturbation  to  the  Instantaneous  steady -state  solution.  Such  a 

solution  would  be  most  useful  In  determining  when  the  heat  transfer  wltti  either 

time -dependent  wall  temperature  or  g  can  be  computed  with  sufficient  accuracy 

(41 

from  quasisteady  relations.  The  author'  '  analysed  the  problem  of  unsteady  free 
convection  In  this  manner.  The  governing  equatlonsi  Equations  (1)  through  (3)* 
are  quite  similar  to  those  used  In  the  free  convection  study  of  reference  (4).  The 
basic  differences  between  the  analysis  of  reference  (4)  and  the  present  analysis 
lie  In  the  boundary  conditions  and,  particularly*  In  the  positions  at  which  they 
must  be  applied.  In  the  analysis  of  the  unsteady  free  convection*  the  boundary 
conditions  are  applied  at  y  a  0  and  y  =  «o .  In  the  present  analysis*  the  boundary 
conditions  must  be  applied  at  y  =  0  and  y  s  6  .  The  film  thickness  6  Is  not 
known  a  priori*  and  it  may  be  obtained  only  as  a  result  of  the  solution.  The 
thickness*  moreover*  varies  along  the  plate  and  with  respect  to  time.  In  order 
to  handle  this  complexity  and  also  to  arrive  at  the  particular  form  of  solutions 
sought*  Equations  (1)  through  (3)  and  their  boundary  conditions  are  first  trans¬ 
formed  in  the  following  manner. 

We  may  define  a  set  of  dimensionless  variables  as 


4 


and 


<1  (*•  *) 


=  ^0  (x»  t),  (x,  t), 


t  )  . . . 

I 


(11) 


where 


Equation  (10)  defines  the  ordinate  q  in  such  a  way  that  the  boundary  conditions  on 
all  the  subsequent  transforxned  equations  may  be  applied  at  17  s  0  and  17  »  1 
only.  It  can  be  seen  from  Equation  (10)  that  the  variable  A  represents  a 
nondimensionalized  liquid  film  thickness  and  is 


2 

A  (^Q*  ...  . . .  )  =»  5(x*  t)  [g{t)/4v  X  J 


(12) 


For  the  time  being,  we  let  the  set  of  variables  |  be  a  set  of  arbitrary 

functions  of  x  and  t. 

The  continuity  equation  (1)  is  satisfied  in  the  usual  manner  by  defining  a 
stream  function  ^  (x,  y,  t)  by  the  equations: 


u 


V 


(13) 


We  now  introduce  a  nondimensional  stream  function  f  and  let  it  be  related 
to  ^  by  the  following  equation: 


^0’  ^1'  ^n’  ' 


^(x.  y.  t) 

— m -  - 

4vx  C(t)A(^Qi  fj*  •••) 


(14) 
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Wo  olso  defino  &  nonditnonsionol  tomperoture  function  oss 

^  ^  y* 

^1*  •••^n'  •••>  •  ~T^tr-  T  • 

w  S 

In  Equations  (10),  (12),  (14),  and  (15),  it  was  assumad  a  priori  that  A  is  a 
function  of  |^^(x,  t)|,  and  f  and  9  are  functions  of  q  and  {f„(x,  t)J  only.  Now 
if  the  transformations 


S  (x,  t )  ^  A  (f  Q»  ^1*  ***^n*  ***^ 

^  (x»  y»  t )  ¥  f(>y,  ^Q,  ,  .  •  •  •  •  • ) 

and 

^(x,  y,  t) . .  ►  6(^1  ^Q»  f j*  •  •  •'  •  •  • ) 

are  carried  out  in  Equations  (2)  through  (9),  and  if  the  resulting  equations  could 
be  made  to  be  functions  of  9  and  |^jj(x*  t)|  only,  then  the  a  priori  assumptions 
wovkld  have  been  proved  to  be  self-consistent  ones  for  the  present  purpose. 

The  definitions  of  nondimensional  stream  function  and  the  variables  yield: 


The  momentum  and  the  energy  equations  are  now  transformed  to  the  following 
equations. 
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Momenttun  equation: 


3f  ^ 


8q2 


8qy 


+  1 


]■»' 


1/A 


1/2 

X 


/AZy  8^f  ®^n  dh  y  8A  ®^n\  4„a2/8£  ^  ®^n 


(17) 
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Energy  equation 


^  li  +  SA^f  51 
^  ^  a»> 


2L^iA^q®l  + 

8q 


(T^-V 

(T^  -  T,) 


-Aij 


(18) 


In  the  above  two  equations,  C  and  (T^  -  T^)'  are  aero  when  the  body  force 
and  the  wall  temperature,  respectively,  are  independent  of  time. 


A  study  of  Equations  (17)  and  (18)  shows  that  all  the  variables  x  and  t  ap¬ 
pearing  explicitly  in  these  equations  would  disappear  completely  when  is 
defined  as: 


n+1 


d-'+l  (T^  .  T,) 


>  lx\  -  1  -  “w 

4-  -(g)  (T,  -  T„l  - jpi 


(19) 


when  T  is  unsteady,  and 
w 


n+1 


1 

.3+2n 


d'>+*C 


dt 


n+1 


(20) 


when  g  is  unsteady.  Here  it  is  considered  that  and  g  are  continuously 
differentiable  with  respect  to  t.  The  fact  that  the  particular  definitions  of 
make  equations  (17)  and  (18)  functions  of  q  and  only,  can  bo  readily  seen 

by  noting  the  relationships; 


8 


■jr 


fof„ 


when  1«  unsteady »  and 


X _ 

vC 


1/2 


TT 


»n+l 


(3  +  2n)^Qf^ 


when  g  Is  unsteady. 

Equations  (17)  and  (18)  now  become  for  the  case  of  unsteady  T^: 
Momentum  equation 

[« ^  -  <i^)'  *  ■]  ■  i:  -  fo  f»)  ^ 

n«0 


1 

8n 


*  ^  2  ^^n+1  "  ^0  dqiV’  ^  ^  9*? 

n*  0  ”  n«  0 


2  8f  ^  ,_..v  ^  8^f 


E  <»**'  fn  -W, 


$  t  t  '-*>  fa  ih] 

n*0  n»0 


2,  « 


(21) 
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Energy  equation 


K  I?  ■  ^ 

*  t  <«„t.  -  fo  41  s  o*‘)  4  il 

n  ’’n  n»0 

nsO 


-  A  £ 


M 

8»? 


E 

n*0 


(22) 


Equation*  (17)  and  (18)  become  for  the  case  of  uneteady  g: 
Momentum  equation 


■^’41:  [4.1- <’^^“'44]^ 

n-O*- 


-  (3+2n) 


t  H  E  <”♦*>  4 

n*0 


a^f 


a»?Hn 


-  2f 


(n+l) 


aA 

»^n 


-  24"  -HE(»+‘>  4 


aq 


af 

9^n 


(23) 
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Energy  equation 


1  B^d 


Pr 


3  £  II  «  ^0  ’'I?  "  I?  ]C  [  ^“+1  '  ^0  ^"1 

n*0  •- 


+  A' 


2  [fntl  -  fo  fnlf^  +  H  S  <"**>  fn  ^ 

n>0  ^  "*  «^0 


H  E  <-*>  t  -  ^  ik 

nsO  n“0 


(24) 


The  boundary  condition*  (4)  through  (8)  are  transformed  for  both  unsteady  -T^ 
and  -g  cases  as: 

At  n  «  0 

e  «  1  (25) 

f  »  0  (26) 

Bl/Bn  .  0  .  (27) 

At  »?  *  1 

0  «  0  (28) 

B^f/Brj^  »  0  .  (29) 
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Boundary  condition  (9)  is  transformed  at  i|  «  1  for  unsteady  case  as: 


<<^.1  -  ^0  ^n)  if  "  ^n  if 

n*  0  n*  0  ^ 

.  ,a2  ^  .  ...  f  »f  1  »« 

*  2A  2^  (n+1)  f„  -  j  .  -  pj  -  — 


nsO 


and  for  unsteady  -g  case  as: 


(30) 


-A"  f„  t  A 


^0  QO 

r^n+l  -  ^0  ^nl  -jf-  +  ^n  if 

n»0  *"  J  n  n 


CO 


t  iA'^S  (ntl)  ^ 

^  ^  n. 

nsO 


J_  ^ 

h®  a>? 


(31) 


Equations  (21)  through  (31)  are  functions  of  9  and  only.  It  is  now 

clear  that  the  initial  assumptions  made  on  A,  f,  and  d  are  entirely  self- 
consistent.  That  is,  A  is  a  function  of  |^n(^>  ^)|  »  nnd  f  and  0  are  functions 
of  Tf  and  |^ri(x»t)|  only. 

In  order  to  solve  Equations  (21)  through  (24)  with  their  boundary  conditions, 
the  functions  f,  $,  and  A  are  first  expanded  into  generalized  Taylor  series 
about  the  steady- state  solutions  as  follows: 

K-i,  (g.  f,.  •••  f„.  ••■)  -  rW  f„W)  t  .  •••]••• 


W’l  *  ■■■  +  'oi'”'  *  "■]■■■ 


(32) 
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0  (n»  f  Qi  ^ 


•  •  • 


•**)  -  H(n)  +[^0^0^*'^  +  +  •**j*-- 


0  00 


(*?)  + 


... j... 


(33) 


and 


Ml 


0* 


f,. -C-)  •••]••• 


♦[^oA„o*-  *fofiV  *••']••• 

When  the  Seriee  (32)  through  (34)  are  substituted  into  Equations  (21)  through 
(24)  and  when  the  coefficients  of  1»  (q,  li,  etc.***  are  collected  in  each 

equation,  an  infinite  set  of  perturbed  equations  results.  For  the  case  of 
unsteady  T^,  the  zeroth-order  and  the  first  two  first-order  perturbed 
equations  are: 


F'"  +  (3FF"  -  2F'^  +  1)  «  0 

^  H"  +  3r^F  H'  *  0 

1^"' +  r^(3FfQ"-  6F‘fQ' +  5F"fp)  «  -  2rAjj(4FF"  -  2F'^  +  1) 

^  00"  +  -  srA^FH'  -  sr^i^w 


I  (35) 


(36) 
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fj'"  r^(3F£j"  -  sr'fj'  +  7r"£j)  -  -  zrAjisrr"  -  2F'^  +  i) 

+  2  {r\  -rAo»?F") 

(3r^  -  4F'  ®  j)  =  -  2yT\QTjH'  +  2/^ 

-  lOrAjFH'  -  TT^fjH' 


For  the  caae  o£  unateady  g,  the  aeroth-order  equatlona  are  the  aame  aa 
Equationa  (35)  and  the  firat  two  £lrat-‘order  perturbed  equatlona  are: 


£  +r2(3Ff."  -  6F’frt'  +  5F"f-)  ar^(2F'  +  qF") 


2rAQ(4FF"  -  2F'^  +  1) 


1^«o"  +^^{3F0(j'  -  2F'0q)  -tSh'  -  8/ApFH'  -  Sr^  f^H' 


fj'"  +r  (3F£j"  -  8F'fj'  +  7F"fj)  *  +  T  f , 


I 


-  2rA|(5FF"  -  2F'^  +  1) 


^  +r^(3F^j'  -  4F'0j)  *  -rAo^H'  + 


-  7r  f|H' 
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Boundary  conditions  (Z5)  through  (29)  for  both  unsteady  and  -g  cases 

become: 

At  r;  «  0 

F  -  £q  »  fj  =  0 
F'  *  fjj'  »  fj*  »  0 


H  «  1,  0Q  »  »  0 


(40) 


At  I?  =  i 


F'*  «  fjj"  *  £j"  “  0 


H  »  ^0  =  ©1  =  0 


(41) 


Also  at  9  ■  1,  boundary  condition  (30)  becomes  for  unsteady  -T  case 


_2  1 

3r^r  ■  -  ^  p  8  ...y  H- 


,  ,  c  (T  -T  ) 

srA^F  +  sr^f.  .  -  ^  ^  ^o’ 


,  c  (T  -T  ) 

ZTAq  +  lOrA^F  +  ^  P  — 


(42) 


and  boundary  condition  (31)  becomes  for  unsteady  -g  case 


15 


sr^F 

-r^  +  srApF  +  sr^ 

^^Q  +  lorAjF  f  7r^£j 


Pr 


c  (T  -T  ) 
p*  ■  w* 


H' 


_1_ 

Pr 


c JT  -T  ) 


c JT  -T  ) 
p  ■  w 


(43) 


Equations  (33)  through  (39)  and  their  boundary  conditions  (40)  through  (43)  show 
that  the  solution  of  each  of  the  perturbed  equations  will  depend  on  the  parameters 
Pr  and  Cp(Tg-TY(^)/h®.  Steady- state  solutions  of  various  condensation  problemsi 
when  they  are  based  on  boundary  layer  equations,  depend  on  the  two  parameters, 
Pr  and  Cp(T,-Tw)/hO,  as  is  seen  in  references  (3),  (5),  (6),  and  (7).  With  the 
particular  transformations  used  in  the  present  analysis,  it  is  seen  that  we  may 
evaluate  the  universal  functions,  F,  Iq,  H,  6 q,  F,  Aq.  etc.,  as  functions  of  the 
same  two  parameters  only,  and  we  may  obtain  the  general  unsteady  solutions  by 
the  use  of  series  (32)  through  (34), 

Equations  (35)  are  the  equations  for  the  steady-state  case,  and  some  of  their 
solutions  can  be  deduced  from  the  results  of  reference  (5)  by  properly  re¬ 
lating  the  variables  of  the  present  analysis  to  those  of  the  reference.  Equations 
(35)  through  (39)  are  integrated  in  the  present  study  by  the  use  of  IBM-7090 
digital  computer  for  several  combinations  of  Pr  and  Cp(Tg-Tv^)/h°.  Two 
typical  results  are  shown  in  Figs.  2  through  5. 


HEAT  TRANSFER 

The  heat  transfer  at  the  wall  may  be  now  obtained  by  evaluating  (Ol/Oy)^ 
from  the  results  of  the  preceding  section.  The  relationship  between  the  instan¬ 
taneous  iocal  heat  transfer  q  and  the  hypothetical  instantaneous  steady-state 
transfer  is  derived  by  the  use  of  Equations  (32)  through  (34)  as: 
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where 


'st 


>  .1/4 

(t;^) 


H  ' 

_ w 

r 


(45) 


The  computed  values  of  ^qw'*  ^0*  «fc.  •  are  tabulated  in  Tables  1 

and  2.  The  coefficients  of  and  of  Equation  (44)  calculated  from  the 
tables  are  presented  in  a  graphical  form  in  Figs.  6  through  9«  These  figures 
and  Equation  (44)  together  enable  us  to  quickly  evaluate  the  ratio  of  unsteady 
local  heat  transfer  to  the  hypothetical  Instantaneous  steady- state  heat  transfer 
when  are  sufficiently  small. 


DISCUSSION 

Figures  6  through  9  show  the  parameters  (^Ow'/H^'  *  Aq//^)  end 
(  ^ivir'/Hw'  '  which  will  be  referred  to  hereafter  as  coefficients  of 

and  respectively,  calculated  from  Tables  1  and  2  for  unsteady  -Tw  and 
-g  cases.  As  can  be  seen  from  Equation  (44),  these  parameters  directly 
represent  the  first-order  effect  of  unsteady  -T^  and  -g  on  the  condensation 
heat  transfer.  Behavior  of  these  parameters,  therefore,  will  be  discussed  in 
the  present  section. 


Unsteady  T^  Case 

Let  us  first  investigate  the  relationships  between  the  coefficients  of 
and  ^1,  and  Prandtl  numbers.  Figures  6  and  7  show  that  the  effect  of  un¬ 
steady  -Tw  on  heat  transfer  becomes  greater  as  the  Prandtl  number  is 
increased.  The  reason  is  rather  obvious.  A  larger  Prandtl  number  implies 
that  either  the  kinematic  viscosity  is  higher  or  the  thermal  diffusivity  is  lower, 
or  it  implies  both.  As  is  seen  from  Equation  (12),  a  higher  v  means,  in 
general,  a  thicker  film.  Both  increasing  film  thickness  and  decreasing  thermal 
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diffuaivity  Increase  the  time  lag  in  thermal  response  of  the  film  to  the  varying 
T^.  Therefore!  the  deviation  of  unsteady  heat  transfer  from  the  instantaneous 
steady- state  value  becomes  greater  as  Prandtl  nuntber  is  increased. 

It  is  seen  in  Figs.  6  and  7  that  ~  ^0^)  positive  whereas 

(^Iw'/H^'  -  ^l/T)  is  negative.  Equation  (19)  shows. -that  both  and  are 
negative  when  and  Tw"  &re  positive.  Therefore,  in  view  of  Equation  (44), 
the  first-order  effect  of  Tw*  on  heat  transfer  is  to  decrease  it  from  the  instan¬ 
taneous  steady-state  value  whereas  that  of  Ty^"  is  to  Increase  it. 

It  is  also  seen  in  the  figures  that  the  coefficient  of  (q  increases  approxi¬ 
mately  with  >^r  whereas  the  magnitude  of  the  coefficient  of  increases 
approximately  with  Pr.  This  means  that  only  the  effect  of  T^'  is  Important 
when  the  Prandtl  number  is  low,  and  the  importance  of  the  effect  of  T^"  in 
relation  to  that  of  Ty^'  increases  quite  rapidly  as  the  Prandtl  number  is 
increased. 

In  Table  1  and  Figs.  6  and  7,  only  the  results  for  Pr  ^  1  are  presented. 
The  reason  for  this  is  that  the  coefficients  of  and  for  the  lower  Prandtl 
numbers  were  found  to  be  sufficiently  small  so  that  the  effect  of  unsteady  T^ 
on  heat  transfer  is  practically  negligible  when  Pr«l. 

The  relationships  between  the  coefficients  of  and  ^2  ,  and  the 
parameter  Cp(Tg-Tyr)/h*^  can  be  also  seen  from  Figs.  6  and  7.  It  is  seen 
that,  for  a  given  Prandtl  number,  the  magnitudes  of  coefficients  of  and 
increase  as  the  parameter  cp(Tg-Tw)/h°  increases.  The  film  is  usually 
thicker  when  Cp(Tg-Ty^)/h^  is  larger  (see  F'b  in  Table  1).  The  time  lag, 
therefore,  in  the  thermal  response  of  the  film  to  the  varying  T^  becomes 
greater  as  Cp(Tg-Ty,)/h*^  increases. 

Equation  (19)  shows  that  the  magnitude  of  increases  with  l/(Tg-Ty,). 
Figures  6  and  7,  on  the  other  hand,  show  that  the  magnitudes  of  coefficients  of 
and  increase  approximately  with  ^Tg-T^  .  Hence,  it  is  seen  from 
Equation  (44)  that  the  first-order  deviation  of  heat  transfer  from  the  instan¬ 
taneous  steady-state  value  increases  approximately  with  l/^Tg-Ty,  for  given 
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value!  of  T^'  and  The  first-order  effect*  therefore*  of  the  unsteady 

on  heat  transfer  Increases  as  (Tg-T^)  Is  decreased  for  given  values  of 

and  T^"*  It  is  interesting  to  note*  however*  that  this  dependence  of  the 
unsteady  effect  on  the  energy-driving  potential  is  much  weaker  than  that  found 
in  reference  (4)  for  the  unsteady  convection  problem. 

The  unsteady  convective  heat  transfer  depended  on  (  only  nnd 


w 


TTL 


ijiy. 


where  T^  is  the  undisturbed  gas  temperature. 

Unsteady  g  Case 

The  effect  of  unsteady  g  on  condensation  heat  transfer  is  of  a  different 
nature  compared  with  the  effect  of  unsteady  -T^  discussed  in  the  preceding 
section.  In  the  latter  case*  the  unsteadiness  is  created  at  one  of  the  boundaries* 
and  the  disturbance  propagated  into  the  film.  The  unsteady  g*  on  the  other 
hand,  effects  the  entire  film  simultaneously  since  g  is  uniform  throughout  the 
flow  field.  Moreover,  the  unsteady  -T^  affects  the  temperature  profile 
directly  and  the  flow  field  indirectly  whereas  the  reverse  is  true  for  the 
unsteady  -g  case. 

There  are  two  separate  factors  which  cause  time  lag  in  the  flow  response 
of  the  film  to  the  varying  g.  They  are  the  inertial  and  the  viscous  drags  of  the 
liquid  in  the  film.  Let  us  investigate  the  behavior  of  these  two  factors  and,  in 
turn,  the  unsteady  behavior  of  the  film  with  respect  to  Prandtl  number  and 
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Cp(Tg-Ty,)/h*^.  The  time  lag  due  to  inertia  is  independent  of  Pr  and 
Cp(Tg-T^)/h'^  eince  it  ia  independent  of  the  total  maaa  of  the  film.  The  time 
lag  due  to  viscous  drag,  on  the  other  hand,  depends  quite  strongly  on  Pr  and 
Cp(T,-T^)/h°.  At  the  higher  Prandtl  numbers,  the  viscous  effect  is  felt  by 
most  of  the  film.  As  the  Prandtl  number  is  decreased,  the  viscous  effect 
becomes  confined  to  the  portion  of  the  film  closest  to  the  wall.  This  phenomenon 
can  be  seen,  for  instance,  by  comparing  the  steady-state  velocity  profiles  given 
in  Figs.  2  and  4  for  Pr  «  10  and  Pr  >  0.01,  respectively.  Thus,  the  time  lag 
in  the  flow  response  of  the  film  to  the  varying  g  at  high  Prandtl  numbers 
(Pr  10)  is  predominantly  due  to  the  viscous  drag,  whereas  the  time  lag  at 
low  Prandtl  numbers  (Pr  0.  1)  is  mainly  due  to  the  inertial  drag.  This  fact 
is  clearly  seen  in  Fig.  8  where  the  coefficients  of  and  are  plotted 
against  Prandtl  number  for  Cp(Tg-T^)/h°  a  0.  1  .  It  is  seen  from  Fig.  8  that 
the  magnitudes  of  coefficients  of  and  increase  as  the  Prandtl  number  is 
increased  when  the  Prandtl  number  is  high.  It  is  because  the  viscous  drag 
continuously  increases  as  the  Prandtl  number  is  increased.  At  the  low  Prandtl 
numbers,  on  the  other  hand,  it  is  seen  that  the  coefficients  of  and 
become  practically  invariant  with  respect  to  the  Prandtl  number  since  the  iner¬ 
tial  drag  is  independent  of  the  Prandtl  number.  The  transition  from  the  regime 
of  dominant  inertial  lag  to  the  regime  of  dominant  viscous  lag  is  seen  from 
Fig.  8  to  take  place  between  Pr  %  0.  5  and  Pr  w  5. 

As  mentioned  earlier,  the  no-interface-shear  boundary  condition.  Equa¬ 
tion  (7),  becomes  less  accurate  as  the  Prandtl  number  is  decreased  below 
about  0.  1.  There  probably  would  be  an  additional  time  lag  due  to  the  interface 
shear,  and  this  effect  could  become  non-negligible  at  the  low  Prandtl  numbers. 
The  magnitudes  of  coefficients  of  and  may,  therefore,  begin  to  increase 
again  as  the  Prandtl  number  is  continuously  lowered  below  about  0.  1. 

Let  us  now  consider  the  effect  of  varying  Cp(Tg-T^)/h°  on  the  coefficients 
of  and  As  previously  noted,  the  film  is  generally  thicker  when 

Cp(Tg-T^)/h°  is  larger  for  a  given  Prandtl  number  (as  was  seen  from/^'s 
given  in  Tables  1  and  2).  When  the  film  is  thick,  a  relatively  small  portion  of 
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the  total  liquid  near  the  wall  is  strongly  affected  by  viscous  drag.  On  the  other 
hand,  when  the  film  is  thin,  practically  the  entire  film  is  affected  strongly  by 
the  viscous  drag.  Therefore,  with  the  aid  of  previous  discussions,  it  is  seen 
that  the  magnitudes  of  coefficients  of  and  sbouid  increase  as  Cp(Tg-T^)/h^ 
is  decreased.  Figure  9,  which  shows  the  variations  of  coefficients  of  and 
^  y  with  respect  to  Cp(T^-T^)/h°,  supports  this  argument  *or  Pr  »  10  . 

It  is  seen  from  Figs.  8  and  9  that  the  coefficient  of  is  negative  whereas 
that  of  is  positive.  At  the  same  time.  Equation  (20)  shows  that  and  are 
positive  when  C  and  C"  are  positive.  Hence,  the  effect  of  C'  on  heat  transfer  is 
a  decrease  from  the  instantaneous  steady- state  value,  whereas  the  effect  of  C" 
is  an  increase. 

Finally,  let  us  briefly  consider  the  reduced-g  experiments  which  are  being 
performed  by  various  people  in  connection  with  space  applications.  Conditions 
of  unsteady  g  usually  prevail  during  reduced-g  experiments.  It  can  be  shown 
from  Equation  (20)  that 


Hence,  the  effect  of  unsteady  g  on  condensation  heat  transfer  becomes  magnified 
quite  rapidly  with  decreasing  g  for  given  values  of  g'  and  g".  It  is  therefore 
important  to  keep  the  unsteady  -g  effect  in  mind  when  the  reduced  -g  experiment 
includes  condensation.  The  effect  is  most  pronounced  when  the  Prandtl  number 
is  high  and,  at  the  same  time,  -  T^)/h°  is  low.  It  is  noted  here  that 

relationships  between  and  8  were  found  to  be  the  same  in  reference 

(4)  for  the  unsteady  convection  as  the  present  relationship  for  the  unsteady 
condensation. 

CONCLUDING  REMARKS 

The  heat  transfer  associated  with  unsteady  laminar  film  condensation  was 
analyzed  for  a  vertical  plate.  Time-dependent  variation  of  either  the  uniform 
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wall  tomperature  or  the  g-£orce  field  was  considered.  From  the  theoretical 
treatment,  the  first-order  deviation  of  heat  transfer  from  the  instantaneous 


steady- state  value  was  obtained  in  terms  of  a  set  of  nondimensional  variables 
I and  parameters  Pr  and  Cp(T^  -  T^)/h°.  From  the  results  obtained,  it  is 
possible  to  determine  when  the  heat  transfer  with  either  time -dependent  wall 
temperature  or  g  can  be  computed  with  sufficient  accuracy  from  quasisteady 


relations. 


The  effect  of  unsteady  on  heat  transfer  was  found  to  increase  as  the 
Prandtl  number  l/(Tg  -  T^)and  Cp/h°  are  increased.  The  unsteady  “T^,  how¬ 
ever,  was  found  to  have  negligible  effect  on  heat  transfer  when  Pr  <<  1. 

Unsteady  -g  was  found  to  affect  the  condensation  heat  transfer  substantially 
for  the  entire  Prandtl  number  range  (0.001  ^  Pr  ^100)  considered.  It  was  found 
that  the  effect  of  unsteady  g  is  insensitive  to  the  Prandtl  number  at  low  Prandtl 
numbers.  The  effect,  however,  was  found  to  steadily  increase  with  the  Prandtl 
number  at  high  Prandtl  numbers.  The  unsteady  -g  effect  on  heat  transfer  was 
also  found  to  Increase  as  Cp(Tg  -  T^)/h°  and  g  are  decreased.  Therefore,  it  is 
important  to  keep  the  unsteady  -g  effect  in  mind  when  performing  a  reduced  -g 
experiment. 

Finally,  even  though  the  present  study  was  concerned  only  with  the  vertical 
plate,  it  can  be  directly  applied  to  a  horizontal  circular  cylinder  by  redefining 
the  variables  q  and  ^  manner  analogous  to  that  done  in  reference  (4)  for 

the  unsteady  free  convection  problem. 
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PLATE 


Fig.  1.  Physical  Model  Considered. 
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Fig  9  Variation  of  Coeffic.ents  of  4  'ind  4j  in  Krjualion  (44) 
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Table  1 


Universal  Functions  for  Heat  Transfer  (Unsteady  T^) 
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H 
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-0. 0098 

0.  0079 
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0.  0092 
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-0.0313 

0.0138 

-1.044 
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KiH 
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-0. 0489 

0.  0204 

-1.081 

-0. 5868 
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0.  2651 

-0. 0102 

0.0374 

-1.005 

-0.4652 

0. 0464 

10 

0.  3140 

-0.0167 

0.  0449 

-1. 010 

-0. 6494 

0.  0896 

0.4597 

-0. 0463 

0.  0695 

-1.044 

-1. 327 

0. 3507 

0.  5341 

-0. 0648 

0.  0821 

-1.079 

-1.  716 

0. 5532 

■CTi 

0.  1490 

-0. 0180 

0.  2104 

-1. 005 

-1.470 

0.  4635 

100 

KSiH 

0.  1766 

-0. 0295 

0.  2529 

-1. 010 

-2.052 

0.8938 

0.2579 

-0. 0806 

0.  3907 

-1. 044 

-4. 1  78 

3.473 

■SH 

0.2992 

-0. 1110 

0.4574 

-1. 079 
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Table  2 

Universal  Functions  for  Heat  Transfer  (Unsteady  g) 
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w 
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0.001 

0.  1 

7.  019 

3.  349 

-0.  9035 

-1. 014 

0.  0005 

-0. 0003 

0.  01 

0.  1 

2.  412 

1.  177 

-0.  3278 

-1.  on 

0.  0017 

-0.  0009 

1. 0 
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0.  5616 

0. 3620 
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0.  0205 

-0.  0139 

0.  05 

C.  2651 

0. 6374 

-1.  1 33 

-1. 005 

0. 0466 

-0.  1148 

0.  1 

0.  3140 

0. 5447 

-0.  6891 
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0.  0651 
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0.  5 

0.  4597 

0.  4065 

-0.  2425 

-1.044 

0.  1340 

-0.  1404 

1.0 

0.  5341 

0.  3792 

-0.  1713 
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0.  1748 

-0.  1592 

1  00 

0.  1 

0.  1 766 

0.  9507 

-3.  828 
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-1.  162 
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